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We probe the effects of strong disorder (2.4 < fepi < 8.6) on superconductivity in thin films 
of niobium titanium nitride and titanium nitride by measuring the microwave electrodynamics in 
coplanar waveguide resonators. We find a gradual evolution of the electromagnetic response with 
disorder, deviating from BCS theory, for both materials. This result can be understood as due to 
changes in the quasiparticle density of states, induced by the short elastic scattering length. The 
observations are consistent with a model using an effective pair breaker, dependent on the level of 
disorder. 



It has been a longstanding paradigm in superconduc- 
tivity, that the properties of the superconducting state 
are not affected by disorder [jQ. Consequently, it has 
become a justified practice to expect a fixed critical tem- 
perature T c for superconducting films, directly related to 
a standard BCS gap and quasiparticle density of states, 
irrespective of their large differences in normal state resis- 
tivity. On the other hand, theoretical evidence has been 
developed, that reveals severe deviations from BCS the- 
ory for materials with a resistivity in the range 100 /xficm 
and higher. In materials with such a large resistivity, the 
elastic scattering length I is of the order of the interatomic 
distance. Therefore, it can be expected that localization 
effects become important. From numerical simulations 
it has become clear that - even for homogeneous disor- 
der - eventually an inhomogeneous superconducting state 
will arise, when the disorder is increased [3J- A short 
elastic scattering length enhances both the Coulomb in- 
teraction between electrons [3], and the interference of 
electrons scattering from impurities. These mesoscopic 
fluctuations were shown to grow when approaching the 
superconductor-to-insulator transition [3], and unavoid- 
ably affect the properties of the superconducting state [5] . 

In early tunneling experiments, it was shown that in 
films of granular aluminum the coherence peak in the 
quasiparticle density of states gets smeared out signifi- 
cantly when the amount of disorder is increased [BJ . More 
recently, it was found that increasing the level of disor- 
der leads to a direct superconductor-to-insulator transi- 
tion for various materials [7j. It was found in thin TiN 
films that, close to this transition, the quasiparticle gap 
shows large spatial fluctuations [5], and by similar tun- 
neling experiments in InO^, evidence was obtained that 
on the insulating side of the superconductor-to-insulator 
transition, there is a peculiar gapped state, presumably 
containing localized Cooper pairs [5]. 

From both experimental and theoretical points of view 
it has thus become clear that superconducting films with 
a large resistivity, are unavoidably prone to deviations 
from BCS theory, even when the disorder leading to the 
large resistivity is homogeneously distributed. Proper 



understanding of these intrinsic deviations is interesting 
from a fundamental point of view, as well as crucial for 
a variety of experiments: the interest in highly resistive 
superconducting films has taken a big leap over the last 
years in the development of sensitive detectors [TO], and 
in experiments involving quantum phase slips |llj . or 
where a superconductor is used in combination with a 
large magnetic field [T^HH] ■ In all these applications, 
a proper understanding of the superconducting state is 
necessary to interpret performance of the devices and the 
results of experiments, especially when the nature of the 
electronic states is an important experimental ingredi- 
ent [TS] . 

Until now, most experimental effort in understanding 
the superconducting state at high levels of disorder has 
been focused on magnetoresistance measurements as a 
function of disorder, yielding information close to T c , or 
on position-dependent scanning tunneling spectroscopy, 
showing the local properties of the material. In this Let- 
ter, we use instead the response of a superconducting res- 
onator at microwave frequencies as a probe for the global 
superconducting state at low temperatures T <C T c , al- 
lowing the convenient study of films with a variety of 
resistivities. We study a series of titanium nitride and 
niobium titanium nitride films with varying normal-state 
resistivity, 120 /xficm < p < 550 /if2cm, still away from 
the superconductor-to-insulator transition. For increas- 
ing disorder, we find increasing deviations from BCS the- 
ory. We show that these deviations are closely linked to 
the short elastic scattering length in the films and com- 
pare them to the available theory for strongly disordered 
superconductors. 

Table [I] gives an overview of the eight films that are 
studied in this Letter. The series of TiN films (A-F) 
was grown using plasma-enhanced atomic layer deposi- 
tion (ALD) [TB]. This process is self-limiting, and there- 
fore in principle allows for growing the TiN film one 
monolayer at a time. In practice, the average growth 
rate of the TiN film is 0.45 A per cycle. The films are 
grown with a thickness ranging from 6 nm to 89 nm. The 
structure of the films is of a polycrystalline nature, as de- 
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TABLE I. Parameters of the films studied in this experiment. 
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termined with atomic force microscopy. The grain size 
b increases with increasing film thickness, and is signifi- 
cantly larger for film F, that was deposited on passivated 
silicon instead of on native silicon oxide. From the fact 
that the other parameters of film F do not differ from 
the trend observed in the other films, we infer that the 
grain boundaries are not limiting the electronic transport 
properties. The NbTiN films G and H were deposited us- 
ing DC magnetron sputtering from a Nbo.7Tio.3 target. 
The resistivity was controlled by varying the pressure in 
the sputtering chamber during deposition |16) . 

All free-electron parameters were determined at a tem- 
perature T = 10 K for the TiN films, and T = 40 K for 
the NbTiN films. Resistance, Hall effect, and upper criti- 
cal field were determined on a Hall bar structure for each 
film. From the measured carrier density and diffusion 
constant, values for the disorder parameter k-pl, where 
is the Fermi wave vector, and the elastic scattering time 
r were determined using free-electron theory [16] . kpl in- 
creases as the resistivity decreases, from 2.4 for the most, 
to 8.6 for the least disordered film. The elastic scattering 
length ranges from 2.4 A to 7.3 A, and is of the order of 
the interatomic distance (a ~ 4.1 A for TiN |17|). 

Coplanar waveguide quarter- wave resonators were pat- 
terned from the films using e-beam lithography and reac- 
tive ion etching. The resonators are capacitively coupled 
to an on-chip feed line, which is wire bonded to coaxial 
connectors. The sample is mounted onto the cold fin- 
ger of a He-3 sorption cooler with a base temperature 
of 310 mK. Microwaves from a vector network analyzer 
are fed to the sample through coaxial cables that are 
attenuated and filtered at 4 K. The amplified forward 
transmission S21 of the feed line is recorded as a func- 
tion of temperature and microwave frequency. From this 
transmission spectrum, the resonance frequency of the 
resonator is determined |16j . 

Figure [T] shows two typical traces of the measured res- 
onance frequency fo (circles) as a function of tempera- 
ture, for the resonators made from films D and A. The 
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FIG. 1. (Color online) Measured resonance frequencies as a 
function of temperature, for films D (top) and A (bottom). 
The dashed curves are calculated resonance frequencies ac- 
cording to BCS theory. The solid curves are a fit using a 
broadened density of states (see text for details). The insets 
show the quasiparticle densities of states used for calculating 
the resonance frequencies. The curves for all films listed in 
Table [I] are presented in the supplementary information [16] . 



resonance frequency decreases with increasing tempera- 
ture, which reflects a weakening of the superconducting 
state. Since our films are in the local limit I -C £, with 
£ the superconducting coherence length, we assume that 
the response to microwaves can be described by a com- 
plex conductivity a = a% — ia% . The resonance frequency 
of the resonators is a function of o~2 and geometric fac- 
tors only |16j . The dashed curves are predictions of this 
temperature dependence using the complex conductivity 
as given by Mattis and Bardeen [TH]. To obtain these 
curves, only the resonance frequency at the lowest mea- 
sured temperature was taken as a scaling parameter [16] . 
It is clear from this comparison, that the measured res- 
onance frequency decays faster with temperature than 
predicted, as was already noticed in Ref. [T5] for a sin- 
gle NbTiN film. Moreover, the deviations from theory 
become stronger for increasing disorder. We attribute 
these deviations to an intrinsic, disorder-induced change 
in the superconducting state. In the remainder of this 
Letter, we will analyze the observed deviations in more 
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detail, and compare them with the available theory for 
strongly disordered superconductors. 

The effective attraction between electrons in a su- 
perconductor is the result of a competition between an 
attractive electron-phonon interaction and a repulsive 
Coulomb interaction [20] ■ It is now known that dis- 
order can change the critical temperature of a super- 
conducting film. This change was first attributed to 
changes in the electron-phonon coupling JT], but later 
it has become clear that disorder has a strong effect on 
the effective Coulomb repulsion. This insight was used 
by Finkelstein [3] to explain the decrease in T c of two- 
dimensional [22] MoGe films with increasing sheet resis- 
tance Rq = p/d. In this framework, the critical temper- 
ature is given by 



T c . = 



\f9~ Vfk + y/gd^g 



(1) 



where g — 2Trh/e 2 R\j is the dimensionless film con- 
ductance, and g c — ln 2 (fi,/fcBT c or)/27r is a critical con- 
ductance, with T c o the critical temperature of the non- 
disordered material. This mechanism, where supercon- 
ductivity is destroyed due to the disappearance of attrac- 
tive interaction, is dubbed the fermionic mechanism, in 
contrast to the bosonic mechanism where phase fluctua- 
tions between weakly connected superconducting islands 
destroy global coherence. We do not include phase fluctu- 
ations in our current analysis, since we expect that these 
will mostly influence the superconducting state close to 
T c , and close to the superconductor-to-insulator transi- 
tion, i.e. far from our experimental situation |23j . 

Equation ([!]) describes the critical temperature of a 
uniform superconductor. The short elastic scattering 
length however will enhance localization, and will lead 
to mesoscopic fluctuations of the pairing potential A — 
A(E,1>) @]. The average value of the pairing potential, 
(A), still follows the BCS temperature dependence. Fluc- 
tuations of A impact the quasiparticle density of states 
however [5] and can be described on a mean field level by 
introducing a parameter rj in the Usadel equation |24j 



iE sin + (A) cos 6 - (A) 77 sin cos 9 = 0, 



(2) 



where E is the energy taken from the Fermi energy, 9 the 
pairing angle, and sin 9 and cos 9 are the quasi-classical, 
disorder-averaged Green's functions. The parameter 77 
is a measure for the spatial correlations in A(r). As- 
suming a 3-dimensional superconductor and long-range 
spatial correlations, its temperature-dependent value is 
given by [5] 



V(T/T C ) 



K(T/T C ) 
4ff(3 - 9c) 



(3) 



where K(T/T C ) is a universal function of temperature, 
with K(0) = 1. 



After solving Eq. ^ , the quasiparticle density of states 
is given by A qp = N Q Re(cos0), with N the density of 
states in the normal metal. For increasing values of 77, 
the coherence peak at E = A is lowered and broadened, 
as qualitatively ilustrated in the insets of Fig. [T] A finite 
density of states appears for energies E < (A). In the 
mean-field description used here, for low enough values 
of 77, the density of states stays gapped. Localized states 
will however create a finite but small density of sub-gap 
states [5] [25]. This so-called Lifshitz tail will not be taken 
into account in the rest of this Letter. 

The imaginary part of the conductivity can be calcu- 
lated using a generalized Mattis-Bardeen equation [26], 
that is valid for a superconductor with arbitrary pairing 
angle 9 describing the extended electronic states: 



<T 2 



f™ nu dEg 2 (E,E>)[l-2f(E>)] 
+ f c °°dEg 2 (E',E)[l-2f(E)}, (4) 



with to the microwave frequency, <7 n the normal state 
conductivity, E' = E + hco, f{E) the quasiparti- 
cle distribution function, and the generalized coher- 
ence factor g 2 (E 1 E / ) = Re[cos 0(£)]lm[cos 9{E')] + 
Im[i sin (£?)] Re [z sin (£')]. For 77 = 0, Eq. Q reduces to 
the standard Mattis-Bardeen expression. Qualitatively, 
a broadened density of states as derived from this the- 
ory, can be used to describe the resonator measurements 
presented here. Quantitatively, the comparison fails how- 
ever. Even for film A, where the broadening effects are 
expected to be largest, only a value of 77 ~ 10 -4 is pre- 
dicted from Eq. ([3|, whereas a value of 77 » 0.1 would be 
needed to describe our measurements. 

To make further progress in analyzing the data, we 
modify the approach taken by Feigel'man and Skvortsov, 
by assuming an arbitrary pair breaking parameter a in 
the Usadel equation ([2]), which then reads 



iE sin + A cos 9 — a sin 9 cos = 0. 



(5) 



This parameter describes a spatially uniform breaking 
of time-reversal symmetry, the origin of which needs to 
be determined. We also assume that, in contrast to the 
model discussed above, this pair breaking mechanism in- 
fluences the magnitude of the pairing amplitude A via 
the self-consistency equation 



A = NnV 



fc B e r 



dJ5Im(sin0)[l-2/(J5)], (6) 



where V is a measure for the effective interaction of the 
undisturbed system, and 0d is the Debye temperature. 
Both V and 9d determine the critical temperature T' c for 
a — > 0. For each value of a we determine the temperature 
where A — > 0, and we identify this with the measured 
critical temperature T c . This description is analogous to 
the description of Abrikosov and Gor'kov for the effect 
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FIG. 2. (Color online) Value of the pair breaking parameter a, 
extracted from fits of the resonance frequency, versus disorder 
parameter 1/kvl, for the series of TiN films (circles) and the 
two NbTiN films (triangles). The line is a linear fit to the 
TiN data. The inset shows the calculated reduction of T c as 
a function of a. 

of magnetic impurities [27], or to the description of a 
superconducting wire carrying a current [2"5| . 

We use this analysis for the measured resonance fre- 
quencies, as shown in Fig. [I] (solid curves). We take again 
the measured resonance frequency at the lowest tempera- 
ture as a scaling parameter, and then we adjust a to find 
the best fit. This procedure implicitly varies the value for 
T cQ for each value of a. The density of states (at T = 0) 
that corresponds to the fitted value of a is shown in the 
inset. We find that the value of a increases with the level 
of disorder 1/kpl, as is shown in Fig. [2] For the series 
of TiN films (circles), a has a linear relation to 1/kpl- 
For the two NbTiN films (triangles), the trend is similar: 
increasing disorder leads to a higher value of a. 

As a consistency check, we analyze the critical temper- 
atures of the TiN films. The pair breaking parameter a 
will reduce the critical temperature, as shown in the inset 
of Fig. [2] The measured critical temperature of films A-F 
decreases with sheet resistance, as shown in Fig. [3] (cir- 
cles). The squares in Fig. [3] show the critical temperature 
when corrected for the effects of a. Increased Coulomb 
interaction also reduces the critical temperature, as de- 
scribed by Eq. 0. Taking both reductions into account, 
using a constant value of r = 1.2 fs (see Table [i]), we find 
a value of T cQ that is constant within 15% (triangles in 
Fig. [3]). This is consistent with the self-limiting character 
of the ALD deposition process, from which we can expect 
that all films are metallurgically identical. 

In our analysis, we have only taken into account mod- 
ifications of extended electronic states. Localized states 
will however be present, and will generate a finite den- 
sity of states in the quasiparticle gap 125] . The effect 
of these localized states will be most pronounced at low 
temperatures, where mostly sub-gap states will be ther- 
mally occupied. At temperatures T < 0.15 T c we observe 
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FIG. 3. (Color online) Measured critical temperature T c as a 
function of sheet resistance Rq (circles), for the series of TiN 
films A-F. Squares show the value of T' c , when the reduction 
due to a is accounted for. Triangles give the value for Tea 
when reduction due to Coulomb interaction [see Eq. 0] is 
also taken into account. 

a slower saturation of the resonance frequency than ex- 
pected from our analysis. We tentatively attribute this to 
the presence of localized states.. A description including 
these states will necessarily go beyond the description of 
the electrodynamics by Nam [26 . At these temperatures, 
the relative change in resonance frequency is of the order 
< 10~ 4 . The small deviations therefore do not affect our 
analysis at higher temperatures. 

Let us now consider the possible origins of the pair 
breaking parameter a. From the fact that the reso- 
nance frequency decreases monotonically with tempera- 
ture, and the fact that the measured deviations are com- 
parable for different substrates and different films, we do 
believe that the role of magnetic moments in the dielec- 
tric surroundings of our resonator is negligible |19j . We 
can also exclude the effect of an uncontrolled level of mag- 
netic impurities. There is a very small concentration of 
iron expected in our films (the T1CI4 precursor gas con- 
tains 0.12 mass-ppm of Fe), but since we expect all TiN 
films to be metallurgically identical, this small fraction 
of magnetic impurity would yield a constant value for a, 
contrary to the trend that is observed. A quantitative ex- 
planation of our results might be sought along the lines of 
disorder-induced spin fluctuations [23||30]- Alternatively, 
the small elastic scattering length I might enhance the ef- 
fect of a minute fraction of magnetic impurities [31) . In 
both cases, the effects on the superconducting state are 
primarily caused by the level of disorder, and are intrin- 
sically present in the materials. 

In conclusion, we have shown that the superfluid part 
of the complex conductivity, 02, as probed with mi- 
crowave resonators, increasingly deviates from conven- 
tional Mattis-Bardeen theory for thin films with normal 
state resistivities exceeding 100 /ificm. We have analyzed 
these results with the conjecture of a disorder-dependent 
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change of the density of states, including a lowering of the 
coherence peak and the occurrence of states for energies 
below the pairing potential A, as suggested by theoreti- 
cal predictions [2 and previous experimental results [5]. 
Using a heuristic approach, inspired by Feigel'man and 
Skvortsov [S], we provide a consistent interpretation of 
the observed electrodynamics of films with these high 
resistivities. Although further details need to be clar- 
ified, these observations underline that disorder has a 
much more pronounced effect on the microscopic elec- 
tronic properties than is often assumed for s-wave super- 
conductors. 
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SUPPLEMENTARY MATERIAL FOR "STRONGLY DISORDERED S-WAVE SUPERCONDUCTORS 

PROBED BY MICROWAVE ELECTRODYNAMICS" 



E. F. C. Driessen, P. C. J. J. Coumou, R. R. Tromp, P. J. de Visser, and T. M. Klapwijk 



FILM DEPOSITION DETAILS 



Films A-F and H were grown on high-resistivity (p > 1 kfkm) Si(100) substrates. Substrates F and H were 
hydrogen passivated by dipping in hydrofluoric acid prior to the deposition. In the case of films A-E, this treatment 
was not performed and there was a thin native layer of silicon oxide present underneath the film. Film G was grown 
on A-plane sapphire. 

The titanium nitride films were grown using plasma-assisted atomic layer deposition (ALD) pQ. In this process, 
TiCU and a plasma of H2 and N2 are stepwise repetitively introduced into an Oxford Instruments reaction chamber 
and react into TiN and gaseous HC1. After each step, the residual gas is removed from the reaction chamber by an 
argon gas purge. During the deposition, the substrate is heated to 400 °C, while the walls of the reaction chamber 
are kept at 80 °C. The average growth rate of 0.45 A was determined with a step profiler on thick films. 

The NbTiN films F and G were deposited using DC magnetron sputtering, at a DC power of 300 W, using an argon 
flow of 100 seem and a nitrogen flow of 4 seem. The pressure in the chamber was varied to change the resistivity of 
the resulting film (4 mTorr for film G, 12 mTorr for film H) 0. 



DETERMINATION OF MATERIAL PARAMETERS 



The thickness d of films B-E is determined with a step profiler on a co-deposited lift-off structure. The thicknesses 
of the other films was estimated from known deposition rates. 

The grain size b was determined from an atomic force micrograph of the surface of the films. To estimate the grain 
size, the diameter of a few representative grains was determined and averaged. 

To determine the electric parameters of the films, we measured sheet resistance Rq, Hall voltage Vn, and (for 
selected films) upper critical field i? C 2,_L f° r the films in Table I of the main article. From these measurements, 
resistivity p, carrier density n, and diffusion constant D were directly determined: 

P = Had, (7) 

(«) 



V R ed 7 



7re \ dT 



T=T C 



with d the film thickness, / and B applied DC current and magnetic field, respectively, and T the temperature. 
To estimate values for the different listed parameters, we use the D rude- Sommerf eld model for the resistivity [3] 

* 

m 

p = —, 10 

ne z T 

with m* the effective mass of the charge carriers, and r their relaxation time. Furthermore, we assume relaxation by 
clastic scattering such that r = 1/vf, with I the elastic scattering length, and «f the Fermi velocity. 
The Fermi wave vector is estimated from the carrier density n, assuming again free electrons, 

k F = ^37r 2 n, (11) 



and the Fermi velocity and wave vector are related through 



«f = — f • 12 

m* 



Combining these equations, a value for / and kpl can be estimated as listed in Table I, from resistivity and carrier 
density only. 
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To estimate the scattering time, the knowledge of I is combined with the diffusion constant that is measured from 
the upper critical held. In a 3-D geometry, it is given by D = vpl/3, yielding 

I 2 

^ = w < 13 > 

The critical temperature was determined from a measurement of the resistive transition. The resistive transition is 
smooth for all films. Therefore, the temperature at which the resistance was reduced by 10% was taken as an upper 
bound of the critical temperature. For most films, the measurement accuracy of T c was 0.1 K due to a resolution 
limit in the thermometer used. 



RESONANCE FREQUENCY OF A COPLANAR WAVEGUIDE SUPERCONDUCTING RESONATOR 

For a thin superconducting him in the dirty limit, with a = o\ — iu^ and (72 > ffi, the sheet inductance at angular 
frequency uj is given by [I] 

L s = ^oAcoth j, (14) 

with A = 1/ 'y ' [Iqloui the magnetic penetration depth, and d the him thickness. 

The resonance frequency fo of an uncoupled quarter-wave coplanar waveguide resonator is determined by the kinetic 
inductance of the superconducting him, and geometric parameters 

fo = ; 1 (15) 

4y(L g + 7 L s )C' 

with I the resonator length, L g the geometric inductance, 7 a geometric parameter converting sheet inductance to 
kinetic inductance, and C the geometric capacitance per unit length. The resonator length and geometric capacitance 
are mere prefactors in this equation, and in our comparison, they are incorporated in a scaling factor at the lowest 
temperature. 



TABLE II. Geometric parameters of the films used in the analysis of our experiments. 



Film 


d 


S 


w 


L s 


7 




[nm] 


H 


M 


[nH/m] 


[mm -1 ] 


A 


6 


10 


2 


311 


288 


B 


11 


10 


2 


311 


273 


C 


22 


10 


2 


311 


256 


D 


45 


10 


2 


311 


238 


E 


89 


10 


2 


311 


220 


F 


55 


3 


2 


437 


482 


G 


300 


3 


2 


437 


391 


H 


50 


3 


2 


437 


486 



The geometric inductance for a coplanar waveguide geometry with central line width S and gap width W, is given 



by 



Mo K(k') 
4 K(k) 



(16) 



where K is the complete elliptic integral of the first kind, k = S/(S + W), k 2 + k' 2 — 1. The geometric parameter 7 
is calculated with [S] 



7 



(S + 2W)< 



32K 2 {k) W{S + W) 



S W 



S Xn \5 S + W I ' S + 2W 



In 



S + 2W W 
8 S + W 



(17) 



with (5 = d/[47rexp(7r)]. Table [TT| lists the values for the different geometric parameters that were used in our analysis. 



The quarterwave resonator is capacitively coupled to a coplanar waveguide transmission line, of which the power 
transmission S21 is measured as a function of frequency. Figure [4] gives a typical resonance curve. The resonant line 
shape is asymmetric due to an impedance mismatch of the transmission line to the 50 coaxial cables. The resonance 
frequency f is extracted from this line shape by fitting the relation [5] 



S21 = 1 - 



Q/Q c 



(18) 



where Q is the loaded quality factor of the resonance, Q c is the complex coupling quality factor, and / is the microwave 
frequency. A small linear background transmission was accommodated in the fit to account for less-than-perfect 
calibration of the cables. The dashed curve in Fig. [4] is a fit of this line shape to the shown resonance. 



-20 




2.815 



2.825 

/ [GHz] 



2.835 



FIG. 4. Resonance curve of a typical resonator. The dashed curve is a fit to the line shape [Eq. ( |18[ )] used to extract the 
resonance frequency fo- 



DATA FROM OTHER FILMS 



Figure [5] shows the measured resonance frequencies for all films presented in the paper. The dashed curves are 
BCS predictions of the temperature dependence, the solid curves are fits using the model described in the paper. We 
note here, that decreasing the geometric inductance Lg would shift the calculated Mattis-Bardeen curves towards the 
measurements. This shift is however much smaller than the observed deviations, even when the geometric inductance 
is ignored. Moreover, this would not account for the observed trend with disorder. 
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FIG. 5. Measured resonance frequencies as a function of temperature, for all films. The dashed curves are calculated resonance 
frequencies according to Mattis-Bardeen theory. The solid curves are a fit using a broadened density of states (see main text 
for details). The insets show the densities of states used for calculating both curves. 



